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I. INTRODUCTION 



If, in addition to the standard model fields, a new pseudoscalar field exists and couples to hyper- 
charge topological number density, it can exponentially amplify hyperelectric and hypermagnetic 
fields in the symmetric phase of the electroweak plasma, while coherently rolling or oscillating. We 
present the equations describing the coupled system of a pseudoscalar field and hypercharge electro- 
magnetic fields in the electroweak plasma at temperatures above the electroweak phase transition, 
discuss approximations to the equations, and their validity. We then solve the approximate equa- 
tions using assorted analytical and numerical methods, and determine the parameters for which 
hypercharge electromagnetic fields can be exponentially amplified. 
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The origin of the cosmological baryon asymmetry remains one of the most fundamental open questions in high 
energy physics, in spite of the effort and attention it has attracted in the last three decades. In 1967 Sakharov noticed 
£SJ \ |jj that three conditions are essential for the creation of a net baryon number in a previously symmetric universe: 
1) baryon non-conservation; 2) C and CP violation; 3) out of equilibrium dynamics . Since then many different 
hypothetical cosmological scenarios in which the three conditions could be fulfilled have been proposed as possible 
scenarios for baryogenesis. 

Among the different scenarios the electroweak (EW) scenario plays a leading role. It is particularly appealing 
because it involves physics that can be experimentally tested in the working colliders and those that will turn on 
during the coming years. Non-perturbative sphaleron processes, at thermal equilibrium at temperatures above the 
EW phase transition, erase any previously generated baryon excess along the B — L = direction. In addition, if 
some asymmetry is generated during the transition it is erased immediately after the phase transition completes by 
t ne same sphaleron processes, if the transition is not strong enough to effectively suppress them Q. The strength of 
the EW phase transition has been extensively studied in the Standard Model (SM) and its popular extensions ||[|, 
'-^ including leading quantum and thermal corrections to the finite temperature effective potential. In the SM the phase 
, transition seems to be second order or even completely absent for those large values of the higgs mass that have 
not been ruled out by LEP II experiment. In addition it became clear that the mechanisms that were considered 
had difficulties to generate enough asymmetry to explain the observed baryon to entropy ratio One of the most 
dramatic conclusions that emerged from these studies was that either new physics beyond the SM able to change the 
character of the phase transition and generate enough asymmetry was relevant at the EW scale, or that new physics 
at much higher energies was responsible for a generation of an asymmetry B-I^O that has survived until today. 

It has been recently noticed || that hypermagnetic (HM) fields could be significant players in the EW scenario 
for baryogenesis. Long range uniform magnetic fields could strengthen the EW phase transition to the point that 
it is strong enough even for the experimentally allowed values of the SM higgs mass. The reason is that only the 
projection of the hyperfields along the massless photon can propagate inside the bubbles of the broken phase, while 
their projection along the massive Z-boson cannot propagate. This well-known effect in conductor-superconductor 
phase transition |Q adds a pressure term to the symmetric phase which can lower the transition temperature. A 
detailed study of this effect in the SM phase transition has been attempted in several recent papers [|| . The results 
are not quite conclusive at the moment, however, this effect could save a baryon asymmetry with B — L = generated 
during the phase transition from erasure by sphaleron transitions in the broken phase and, therefore, may fix one of 
the two main SM dissabilities discussed above. 

A subtle effect of hyperelectromagnetic (HEM) fields in the EW scenario may also solve the other basic problem 
for EW baryogenesis, the amount of asymmetry that can be generated. Giovannini and Shaposhnikov have shown j(| 
that the topological Chern-Simons (CS) number stored in the HEM fields just before the phase transition is converted 
into a fermionic asymmetry along the direction B — L = 0. 

We have shown ||, that an extra axion-like pseudoscalar field coupled to the hypercharge topological number 
density can amplify HM fields in the unbroken phase of the EW plasma, while coherently rolling or oscillating around 
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the minimum of its potential. This mechanism is capable of generating a net CS number that can survive until the 
transition and then be converted in a baryonic asymmetry in sufficient amount to explain the observed baryon to 
entropy ratio. 

Pseudoscalar fields with the proposed axion-like coupling appear in several possible extensions of the Standard 
Model. They typically have only perturbative derivative interactions and therefore vanishing potential at high tem- 
peratures, and acquire a potential at lower temperatures through non-perturbative interactions. Their potentials take 
the generic form VqV(4>/ /), where V is a bounded periodic function characterized by two mass scales: a large /, which 
could be as high as the Planck scale, and a much smaller mass m ~ Vq /f, which could be as low as a fraction of an 
eV, or as high as 10 12 GeV. A particularly interesting (pseudo)scalar mass range is the TeV range, expected to appear 
if potential generation is associated with supersymmctry breaking. Scalars with axion-like coupling to hypercharge 



electromagnetic fields were previously considered in |l0| , 
Amplification of ordinary electromagnetic (EM) fields 



jy a scalar field was discussed in |110j . In this paper we discuss 
in detail the dynamics of the pseudoscalar field and how it drives HEM fields in the EW plasma at temperatures above 
the phase transition. In section II. we present the basic equations and several useful approximations to them. An 
analytical description of the amplification of HEM fields by this mechanism in the different regions of the parameter 
space is presented in section III. In section IV we compare this analytical approach with numerical results. Our 
conclusions are summarized in section V. 



II. ELECTRODYNAMICS DRIVEN BY AN AXION-LIKE PSEUDOSCALAR. 



In this section we discuss hypercharge electrodynamics in the unbroken phase of the EW plasma coupled to a 
cosmological pseudoscalar field. We assume that the universe is homogeneous and isotropic, and can be described by 
a conformally flat metric ds 2 = a 2 (rf)(drj 2 — dx\ — dx\ — dx 2 ) where a(rj) is the scale factor of the universe, and r\ is 
conformal time related to cosmic time t as a(rf)drj — dt. 

In addition to the SM fields we consider an extra pseudoscalar field <j) with coupling to the U(l)y hypercharge field 
strength jifiYY. The scalar field rolls or oscillates coherently around the minimum of its potential, at temperatures 
T > 100 GeV just above the EW phase transition, and drives the dynamics of the HEM fields. The coupling constant 
A (which we will take as positive) has units of mass -1 A ~ 1/M. For QCD axions, the two scales M and /, which is 
the typical scale of variation for the scalar fields, are similar M ~ /, but in general, it is not always the case. As we 
are not considering here any particular extension of the Standard Model we take M to be a free parameter, and in 
particular allow / > M. 

We also assume that the universe is radiation dominated already at some early time ij = 0, at T > lOOGeV, 
before the scalar field becomes relevant. Once the scalar field starts its coherent motion it can become the dominant 
source of energy in the universe. We will assume that after a short period of time the scalar field decays, perhaps to 
(hypcr)photons or other particles, so the radiation dominated character of the universe is not significatively altered. 

The lagrangian density describing HEM fields, coupled to the heavy pseudoscalar in the resistive approximation 
|12j of the highly conducting EW plasma is 

C = ^—g Qv^W - V(<f>) - -Y^Y"" - W - -^Y^ - ^e ijk Y»Y k , (1) 

where g = de^g^) is the determinant of the metric tensor, V M denotes the covariant derivative, and J M is the ohmic 
current. 

The last term in (Jlj) takes into account the possibility that a fermionic chemical potential survives in the unbroken 
phase of the EW plasma . A finite right electron density could survive if the higgs mediated electron chirality flip 
interaction is sufficiently suppressed, due to the smallness of the electron Yukawa coupling, so that it cannot reach 
thermal equilibrium (^JT^]. The coefficient [i in |l]) is related to the right electron chemical potential fiR by p, = j^z^r, 
where g' is the hypercharge gauge coupling constant. The latin indexes k) run over the spatial degrees of freedom 
only. 

The equations of motion that follow from lagrangian (Q) have been partially derived by several authors [|||l0|,|l5| . 
The complete set of equations of motion includes an equation for the pseudoscalar, 

v^ + ^ = -^y- ( 2) 

an equation for HEM fields, 
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(3) 



and the Bianchi identity, 

v Al r A " y = o. (4) 

We have introduced the four-vector S a = (//, 0, 0, 0) in order to keep the notation clear. 

A plasma is a highly conducting ionized medium, so that individual (hyper)electrical charges are exponentially 
shielded ~ e~ r / TD over distances longer than the Debye radius rr>, which in the EW plasma is of the order of the 
inverse of the temperature of the plasma T, r^ 1 ~ (10 — 100)T Therefore, we will assume that the plasma is 
electrically neutral over length scales larger than rjj. Since we are interested in the coherent motion of the time- 
dependent scalar field we will assume, in addition, that the spatial derivatives of <f> are negligible compared to the 
other terms in the equations. 

Maxwell's equations for HEM fields then become 

(i) V-B = 
(Hi) V • E = 

(iv) ^=VxB + ^ma(v)B - X^B - J, (5) 
or\ 7r z ar\ 

where V now represents the usual three-space gradient (for comoving coordinates). The current J is given by Ohm's 
law in a plasma whose bulk motion is described by a non-relativistic velocity field v, 

J = a(E + v x B). (6) 

In equations (|^,^) we have introduced rescaled electric and magnetic fields E = a 2 (rj)£, B = a 2 (r])B, and physical 
conductivity a = a(rj)a c . The fields £ , B are the flat space HEM fields defined through 
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Equation (^) for the scalar <fi becomes, after expanding the covariant derivative, 

2aH^f- + a 1 - ' Y' =\a 2 E- B. (8) 



drj 2 dq 



The Hubble parameter H = ^ is related to the temperature T, H = -p-, where M Q — 1 ^y= > and is the 
number of effective degrees of freedom in the thermal bath of the plasma. At the EW scale = 106.75 for the 
SM degrees of freedom in the unbroken phase. For us, although the difference is not very relevant in our analysis, 
.9* = 107.75 if we consider the extra bosonic degree of freedom. 

In the absence of the extra pseudoscalar field, Maxwell's equations in a plasma can be solved in the magnetohydro- 
dynamics approach (MHD) of infinite conductivity a — *■ oo ]lj], 

E^-vxB. (10) 

In this approximation the electric E and magnetic fields B are perpendicular to each other, so that E ■ B = 0. We 
can therefore solve eq.(|8|) with a vanishing r.h.s., as a first approximation, and substitute the resulting (f>(rj) into 
eq.^.w) as a background. At the end of the section we give an estimate of the error coming from neglecting the HEM 
backreaction in the scalar field equation. 
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The explicitly time dependent background spoils, however, the validity of MHD approach to Maxwell's equations. 
We therefore need a generalized ansatz to describe the solutions to these equations. Our guess is that the HE field is 
proportional to the time derivative of the HM field, 

-± a dB -±„ , . 

E = -— + ve, 11 

a or] 

where a is a dimensionless coefficient to be determined. We add the term V0 for generality and consistency. In the 
choice of this ansatz we have been guided by two requirements: 

1) when the scalar field decouples, the new ansatz should be the leading correction in the small parameter — to the 
MHD solution; 

2) the HE field E depends linearly on the HM field B. 

Equations (^.i) and can be understood as constraints on the initial conditions, which are conserved in time 

when E and B evolve according to the other two equations. Inserting ( p"l| ) into ^.ii) we obtain 

VxB = --B (12) 
a 

V x E = --E +-V0. (13) 
a a 

The last relation is consistent with (^.m), only if we require, in addition, V 2 6> = 0. Now we have the same constraint 
on the spatial dependence of both B and E fields. This constraint can be most easily solved by assuming a solution 



with factorized time dependence. The spatial modes that satisfy eq.(12) can be labeled by their wave number k and 
a sign ± for the two possible helicities, 



with 



B% = e- ik %±P±{ri), (14) 



bz ± = b±(e 1 ±ie 2 ). (15) 



Here ei, e.2 are unit vectors in the plane perpendicular to k such that (ei, £2, k) is a right-handed system. When we 
compute the curl of these spatial modes we find that the proportionality constant a in our guess ( pi] ) , for each given 
mode, should be 

a = ±T- (16) 
k 

Once we have identified the magnetic modes we can, in a straightforward way, obtain the corresponding electric 
modes from cq. (|lT|), 

% = e ~ iSx V e fM+™, (17) 



with 



k ^e±(r 1 )=±b^ ± . (18) 



drj 



When we substitute expressions 

(HJl^) into eq-d-iu) 

we obtain a second order equation for the time dependent 



factors (77) , 

Equation ( |l9| ) is the main result of this section. The rest of the paper is devoted to finding useful approximations to 
it, analyzing their solutions, and to numerical studies of its solutions. 

The term ± 177 k fit (rj) is the so-called dynamo-term. Such a term is not unusual in plasma physics. It has been 
claimed |l7j that such a term could give rise to the long-range homogeneous magnetic fields in astrophysical systems. 
A similar effect is provided by the /i-term H . The coefficient 7 is defined in the approximation 
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V x (vx B) w 7 (V x B), (20) 

according to the procedure outlined in p7| . If the bulk motion of the plasma is random and has zero mean velocity 
< v >— 0, then it is possible to average over the possible velocity fields, assuming that the correlation length of the 
magnetic field is much larger than the correlation length of the velocity field. The procedure is equivalent to averaging 
over scales and times exceeding the characteristic correlation scale and time To of the velocity field. On the other 
hand, if the characteristic correlation scale and time of the velocity field are much larger than those of the magnetic 
field we can assume v = 0. The coefficient 7 is given by 

7 = -|(fVxS), (21) 



and so, vanishes in the absence of vorticity in the plasma bulk motion. From eq. (|20|) it is possible to conclude that 

vxBkjB + Wuj (22) 

for a certain given scalar function uj. 

The consistency of the analysis that we have presented imposes 

— — + aVe + cjVlu = 0. (23) 
or\ 

This constraint can be obtained if we subtract (|^.») from (|B|.w) and then insert our particular solution. In the limit 
of infinite conductivity a — > 00 and in the absence of vorticity (7 = 0) we have 

V0 = -Vw = -(vx B) + jB = -vx B. (24) 

In the particular case when the extra pseudoscalar is not present our ansatz describes exact modes for electrodynamics 
in a conducting medium. 

When the extra pseudoscalar is present, an exact analysis should take into account the backreaction of the electro- 
magnetic modes on the scalar equation (H). We can estimate the error in dropping this backreaction. 
The concrete scalar potential that we use for the analysis is 

V(4>) = V 4 V(q>/f) = m 2 f 2 V(^/f) » im 2 ^ 2 , (25) 

for cj) < /. For values of m in the TeV range and T Z lOOGeV, the mass term in ( dominant 
over the cosmic friction term 2aH^. If we drop the electromagnetic backreaction, the scalar field equation is that of 
an harmonic oscillator whose period^ 2n/m, is much shorter than the characteristic expansion time of the universe 
at that time, jj ~ Mpi anc k/T 2 . Consequently, we can fix a(rj) = 1 for the period of time that scalar oscillations last. 
The general solution for the scalar field in this case is 

4>(r]) = Acos(mr] + po) ~ fcos(mrj + po), (26) 

where po is to be fixed as an initial condition. 
The backreaction term is negligible if 



\\E-B\<&. (27) 



In a previous paper Q we have found, under very general assumptions and using our particular ansatz to describe 
the electrodynamics driven by an oscillatory scalar background, that after averaging over magnetic and scalar fields 
initial conditions, (E ■ B) can be expressed as follows 

^■B) = ^^-(BLJ-KJ) (28) 
where B^ is a particular mode that is maximally amplified. Below we will show that k max is approximately 
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If we integrate over a scalar oscillation we obtain average values over time. Backreaction is negligible if 



Multiplying both sides of the equation by fm and dividing by \{B^ 2 + B k 2 )|, we obtain 



« 77^Z ; . (31) 



2k max /m\ { BlJ + B- m J)\ \(Bt m J+B kmax „ 

The right hand side of the inequality is less than unity if the energy density in the amplified magnetic modes exceeds 
that stored in the scalar oscillations. So the condition can be estimated as 

\(B+ 2 ~Bl 2 )\ 2k max /m , s 
J mal < °f w 1, (32) 

\( B LJ + B ^J)\ f 

according to (p9|). 

The asymmetry parameter 

(B+ 2 -B7 2 ) 

= \^ ^± (33) 

can vary from - no asymmetry - to ±1 - total asymmetry. As we will show below (see also ||), in the general case 
of an oscillating scalar field both modes (±) are roughly amplified by the same exponential factor, so that jb remains 
approximately equal to its initial value. There we also showed that this asymmetry parameter may be related to the 
amount of baryonic asymmetry generated at the EW phase transition. We concluded that only when this parameter 
is significatively small the generated asymmetry is consistent with the experimental ns/s ~ 10~ 10 . This justifies why 
we may neglect the electromagnetic backreaction in the scalar equation. 



III. MAGNETIC FIELD AMPLIFICATION: ANALYTIC DESCRIPTION 



Before trying a detailed description of the general solution of eq.(19), it is useful to consider the simpler case of a 
constant ^ . The solutions are then simply linear superposition of two exponentials 



where lo 1 2 are the two roots of the quadratic equation, 



which are 



k 2 ± A— k ± ajk =p ^ n R ka(rf) 

CtT] 7T Z 



(34) 



(35) 



'1,2 



-a± 



4 [k 2 ± \-?-k ± 0-7/c =f ^-zV>Rka(rj) 
dr] 7H 



(36) 



Unless a particular choice for the initial conditions is made, the solution is largely dominated by the exponential 
which corresponds to the larger of the two eigenvalues. This is an exponentially growing function if uj\ is positive. 
The linear terms in the wave number k can exponentially amplify one of the two helicity modes for a limited wave 
number values. For larger wave numbers the quadratic term k dominates over the linear terms and the magnetic 
modes are oscillating and/or exponentially damped, depending on the overall sign of the expression under the square 
root. 

The right electron chemical potential ]Xr can be computed j^] through the expression 



2 2 783 201 15 
m = 45" 9 * [ W 6R ~ + 22 {h + 5s)]T > 



(37) 
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where Sr = nn/s (s is the entropy density) is the right electron asymmetry. For generality we have included the 
possibility of finite asymmetries <5, in the three conserved charged Lj — ^. The index i = 1, 2, 3 runs over the three 
fermion families. If we assume that these asymmetries can be of the order of the observed baryonic asymmetry 
Ub/s ~ 10" 10 O] we get that the linear term in k, ^rP-Ra(r]), in eq. ( |36| ) is of the order 10~ 10 T, where T is the 
temperature of the plasma. In the model that we are exploring the scalar velocity term ~ A/m > m is much 
larger than the chemical potential contribution. We will then fix for simplicity pn = in the following. 

We will assume that the velocity field v is parity invariant, so that no vorticity at the length scales A: -1 we are 
considering is present. We then fix 7 = in addition to p = 0. Therefore, amplification can occur for one of the two 
helicity modes if 



dr) 



> k. (38) 



To obtain significant amplification, coherent scalar field velocities ^ over a duration are necessary, larger velocities 
leading to larger amplification. 

It is interesting to remark at this point that, in absence of the scalar field, and in the limit a — > 00, the magnetic 
modes (both helicities) diffuse, 

B k (v) ~ (39) 
as MHD predicts §. Therefore, (see eq. (0) and (H)) 

-. 1 dB - k -> 

E(rf) = -vx B ±-— « -vx Bt -B, (40) 
k OTj a 

that confirms, as we had anticipated, that in the appropriate limit our exact modes describe the leading correction in 
1/cr to the MHD approach @. 

In the general case of an oscillating scalar field (^6|), the field's velocity changes sign periodically and both modes 
can be amplified, each during a different part of the cycle. Let us first redefine eq. ( p^ ) in terms of the dimensionless 
parameter u — mr], 

d 2 ^ , a dfc* ( ( k\ 2 , A . . .(k 



+ ^if + [{™) ±Aam ( u + , )^-jJ^W = 0, (41) 

where we have introduced the dimensionless parameter A = XA. We can expect, on the basis of the example described 
at the beginning of the section, that each of the two modes will grow exponentially during parts of the scalar oscillation 
when 



± Asin(u + p ) < 0, (42) 



fe 
m 

respectively, and oscillate or be damped during the other part of the cycle. Equation (^2|) imposes an upper bound 
on the wave number spectrum that can get amplified, since —1 < sin(u) < +1, we need that 

k < Am. (43) 

Net amplification results when growing overcomes damping during a cycle, total amplification is then exponential 
in the number of cycles. In order to estimate in which cases this will happen, let us assume that k is small enough, 
k/m < A. Then, according to ([i"2"|), each mode gets amplified during half a cycle and damped during the other half. 
If, in addition, 4 |±Asm(u + /°)(m)l < ^(m) > ^ ne re l evan t eigenvalue lo\ can be approximated by a linear expression 

wf w - (±Asin(u + p)-^j , (44) 

so that the growth during half a cycle is canceled by the damping when the sine changes its sign, and therefore no net 
amplification results. But, if the two terms in the last inequality are at least comparable, the linear approximation is 
not valid; a quadratic correction, that does not change its sign with the sine function, is relevant and net amplification 
results. This gives us a lower bound on the amplified wave number spectrum, 
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A 



m I \m 



(-Y- (45) 



We conclude that amplification occurs for a limited range in the wave number spectrum. Maximum amplification 
occurs for that value of k that maximaizes the expression (^) ± A f^J sin(u), namely, 



|m. (46) 



To be more precise, this estimate is an upper bound on k max . We have not taken into account that the larger k is, 
the smaller the duration in which amplification occurs (see condition (H3)). Nevertheless, numerical analysis confirms 



that fl46D is quite accurate. (For example, see Fig. 



The conductivity of the plasma a, as well as the comoving wave number k are both proportional to the temperature 



of the plasma T, a — oqT ~ 10T |16|] and k = k$T. The character of the solutions of eq.(41), for a given comoving 
wave number, depends on the dimensionless parameter T/m. In Fig. |3|we show that for a given ko there is a limited 
range of temperatures where the oscillating scalar field drives net magnetic amplification, as can be seen from eqs. 
© and ©. 

The situation is somewhat different when m <C T. In this case eq. (^lj) can be approximated by a first order 
equation 

which can be solved exactly, 



ml du \ m 



P±(t,) = 0$ T K<Piv) ~ 0(0)) (48) 

In this extreme limit, the friction term in the scalar equation (0) is comparable to the mass term. It means that it 
would take a time interval of the order of the characteristic cosmic time expansion for the scalar velocity to change 
its value significatively, so for our purposes, we may simply consider d(j)/dr) ~ constant. The scalar field does not 
oscillate during the relevant time scale, instead we say the scalar field rolls. When the scalar field rolls only one of 
the two helicity modes gets amplified. The hclicity mode that is amplified, is determined by the sign of <j>(rj) — 0(0). 
The amplification factor A (k, n) = fifc (n) / fjfi- (0) due to this mechanism is maximal for the wave number k, that at 
a given time 77, maximizes the exponent in eq. (J4q) : 

k max V = ^AA0, (49) 



A ± (k maxlV ) = £. (50) 

i(0) I . Looking at 

\Acf> > 10 s . So we have 



Here A<p(n) = \<f>(rf) — 4>(0)\. Looking at 77 ~ r\EW we obtain ~ 10 16 , and therefore to obtain amplification 



^max,EW , „8 ir ,-16 ° EW 



. 10 8 • 10- 16 ■ ^ » 10- 7 . (51) 
Tew Tew 

A value of AA<^> > 10 8 is not unnatural, for example, such a value is obtained if the typical scale for scalar field 
motion is the Planck scale, as happens in many models of supergravity, and A < \/\{) w GeV . 

In the rolling case the discussion about the role of the electromagnetic backreaction on the scalar equation (||) is 
somewhat different from the oscillating case discussed at the end of the previous section. As we have just shown only 
one of the two helicity modes is amplified when the scalar field rolls, while the other is damped. The backreaction 
term \a 2 E ■ B can be expressed in terms of the solution ( p0[ ) 

where k max is the mode specified in 
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This term is to be compared to any of the dominant terms in the scalar equation (m. In this case we choose for 
convenience the friction term 2aH§&. The backreaction is negligible while 



A , d 



If we integrate over time then the condition can be written 

4^k max HA(j) 



< 1. (54) 



When bound (|54| ) is saturated the electromagnetic backreaction in the scalar field equation becomes relevant and 
changes the free motion of the scalar field. 

A complete discussion of the end of the coherent oscillations or rolling by the scalar field is beyond the scope of 
this paper. However, we would like to make two comments in this regard, 

• we do know that once the oscillations or rolling stop, the fields are no longer amplified and obey a diffusion 
equation (^9|). Modes with wave number below the diffusion limit k < k a T ~ 10~ 8 T, where = F 
remain almost constant until the EW transition, their amplitude goes down as T 2 , and energy density as T 4 , 
maintaining a constant ratio with the environment radiation. Modes with k /T > k a decay quickly, washing out 
the results of amplification. We have seen that the range of amplified momenta for oscillating fields is not too 
different than T, therefore scalar field oscillations have to occur just before, or during the EW transition, for the 
mechanism we are discussing in this paper to be relevant for EW baryogenesis. In that case, the amplified fields 
do not have enough time to be damped by diffusion. If the field is rolling, momenta k <ti T can be amplified 
and then frozen in the plasma until the phase transition, and therefore the rolling can end sometime before the 
transition. 

• Some clues about how oscillations or rolling may eventually end can be obtained from the estimated HEM 
backreaction term in the scalar field equation. For an oscillating field, we have seen that the backreaction term 
remains negligible throughout the evolution. In that case, another type of effect or interactions have to be 
considered as leading eventually to the end of coherent motion. If the field is rolling we have seen that the 
backreaction term can become significant, and therefore it may well be that in this case the coherent decay of 
the scalar field into HEM fields is the cause for the end of coherent motion. 



IV. NUMERICAL ANALYSIS. 

We have studied the numerical solutions of eq.(|4l]) in different regions of the parameter space. We find results 
that are in very good agreement with the previous qualitative discussion: amplification occurs for a limited range of 
Fourier modes, peaked around k/m ~ \h, (see eq. ([h3]) and Fig. ||). The modes of the EM fields in the spectrum 
range that is amplified are oscillating with (sometimes complicated) time dependence and an exponentially growing 
amplitude, 

flk fa) = e r * SFpfa). (55) 

The coefficient Y gives the amplification rate per oscillation of the corresponding mode. The function P(rj) is periodic 
over a period Arj — lix jm. 

In Fig. ^ we show the time dependence of three representative magnetic modes for specific values of parameters: 
(1) k outside the amplification region; (2) low k values; (3) k values around maximum amplification. 
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Magnetic mode with large k 
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Magnetic mode with small k 

Amplification rate per oscillation = 2.1 
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Magnetic mode with medium k 

Amplification Rate per Oscillation = 27.6 




FIG. 1. Amplification of EM fields. The function /?, exponentially scaled, is shown for the parameters A = 55, m = 6TeV 
and a = 40T at T = ITeV, for three different representative wave numbers. 

This mechanism is a very efficient amplifier of EM fields. For example, to obtain an amplification of 10 12 for A = 55, 
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k/m = 20, m = 6TeV and a = 40T, for oscillations occurring at a temperature of ITeV, we need just one cycle! 
Other examples: for A = 18, k/m = 2, m = 6TeV and a = 70T and oscillations occurring at T = 120GeV, in one 
cycle magnetic fields are amplified by a factor 10 4 ; for A = 50, k/m — 10, m = 6TeV and a = 40T at T = ^OGeT^, 
the amplification factor is 10 17 . 

For the range of parameters in which fields are amplified, the amount of amplification per cycle for each of the 
two modes is very well approximated by the same constant T(k/m, A, a). A good approximate estimate for the 
average amplification after N cycles is therefore (fc, 77) = A/' ± fee Arr , where N^k represents the transient influence 
of the initial conditions of EM and scalar fields. 

In Figures || and we have shown an example of the claimed dependence of T on wave number and temperature, 

Amplification Rate Per Oscillation: 

Dependence on Wave Number 
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FIG. 2. Amplification per cycle of EM fields, as a function of wave number. 
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Amplification Rate per Oscillation: 

Dependence on Temperature 
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FIG. 3. Amplification per cycle of EM fields, as a function of temperature. 

In Fig. ^] we represent the amplification rate per oscillation as a function of the amplitude of the scalar oscillations 
A. We show the range of variation with the conductivity a for range a = 10 — 70T. A notable feature of this picture 
are the dips for certain values of A in both graphs. We believe that these are specific values, for which through the 
coupling of higher Fourier modes of the periodic function P{r\) in (^5|) to the scalar oscillation the leading exponential 
factor is canceled. But we do not have a clear understanding of this phenomenon. 
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Amplification Rate per Oscillation: 



Dependence on Scalar Field Amplitude 
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FIG. 4. Amplification per cycle of EM fields, as a function of scalar field amplitude. 

V. CONCLUSIONS 

A pseudoscalar field coupled to the hypercharge topological density can exponentially amplify HEM fields and 
develop a net Chern-Simons number in the symmetric phase of the EW plasma while it rolls or oscillates around the 
minimum of its potential. This mechanism could drastically change the electroweak scenario for baryogenesis and 
perhaps fix its two main dissabilities: the amount of asymmetry generated by electroweak processes and the character 
of the phase transition. 

We have studied this mechanism for (pseudo)scalar masses in the TeV range, that could naturally appear if the 
scalar field is associated to supersymmetry breaking. In that case the coherent oscillations have to occur just before 
or during the phase transition in order to avoid the fast diffusion in the plasma of amplified magnetic modes once 
the scalar coherent motion terminates. In such a case the amplification spectrum is sharply peaked around the wave 
number k max ~ ^Am, total amplification is exponential in the number of cycles, and large amplification of magnetic 
modes, even 10 12 or larger, can happen just after a few scalar oscillations depending on the particular values of the 
parameters of the model. 

If the scalar field rolls instead of oscillating, the mechanism would be relevant for EW baryogenesis even if it takes 
place at higher temperatures before the phase transition. In this case, modes with wave number much smaller than 
the temperature of the plasma k <C T are maximally amplified. Once the scalar rolling terminates, the amplified 
magnetic modes remain frozen in the plasma and do not diffuse. 

In a previous paper we have concluded that our mechanism would be able to generate enough asymmetry to explain 
the baryon number density to entropy ratio observed in the universe. 
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